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1. Introduction 

Utilizing Feshbach resonance physics, recent experiments have produced an atomic Bose- 
Einstein condensate (BEC) coherently coupled to molecules in high vibrationally excited 
bound states [l]. The coupling of atomic and molecular states was achieved through 
application of a pulsed external magnetic field and it has sparked much interest j2] since 
it ultimately (although possibly not yet) should lead to the creation of a molecular 
condensate. The problem of creating a molecular BEC has been one of the focus areas 
of ultracold physics research for several years now [3^ . Two different techniques, namely 
photoassociation and a magnetic field ramp near a Feshbach resonance, have been used in 
attempts to transform an atomic condensate into a molecular one. The use of Feshbach 
resonances to control the atom-atom scattering length, and other properties, has been 
previously demonstrated experimentally [H E] . 

It was predicted theoretically [Hj that magnetic field pulses would drive a significant 
part of an atomic BEC into a molecular one. These predictions were based on a mean 
field theory approach, of the same type that has proven very successful in describing 
many properties of the alkali atom BECs produced experimentally to date. However, two 
sets of experiments performed at JILA have shown puzzling results that did not match 
the original theoretical predictions. In one experiment, a single magnetic pulse close to 
the Feshbach resonance probed the strongly interacting atomic dynamics, [Zj while in a 
second experiment double pulses generated interference patterns between the different 
states populated. In both experiments, the atoms, part of the initial BEC cloud, were 
observed to end up in one of the following three components: a remnant BEC cloud, 
a burst of hot atoms, and a missing (undetected) component. This outcome differed 
from the theoretical predictions which only accounted directly for two components. 
The questions regarding the nature of the three components observed experimentally, 
including the specific issue of whether a molecular BEC was created, were addressed by 
two independent theoretical papers [Hlin]- Both of these papers accounted for many body 
effects by employing field theory beyond the mean field (i.e., Gross-Pitaevskii equation) 
level. These two studies gave similar answers, in identifying the observed atom bursts 
as hot, non-condensed atoms, and the missing component as a molecular condensate. 
In addition to Refs. |Hj and [H], three other theoretical studies involving many-body 
approaches have also addressed the interpretation of the Donley et al. experiment 

[miiniiE]. 

In this paper we propose an alternative model that does not rely on field theory. 
Our model considers only two-body dynamics and uses a very simple scaling procedure 
to apply our results to the many-body system studied experimentally. The observed 
oscillatory behavior can then be viewed as a simple example of quantum beats of the type 
that arise whenever indistinguishable quantum mechanical pathways associated with two 
or more impulsively-excited stationary states interfere coherentlv. |13j Specifically, our 
treatment considers two atoms confined by a spherically symmetric harmonic oscillator 
potential, which interact through a contact potential. The stationary states of such a 
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system have been described by Busch et al, [T3| who proposed it as a model that can be 
used in the context of ultracold colhsions occuring in trapped gases. The properties of 
this treatment were further investigated in the context of the condensed Bose gases by 
Juhenne and coworkers [Tallin] and by Blume and Greene [T7j. These papers investigated 
the strong interactions that occur near a Feshbach resonance, and employed an energy- 
dependent scattering length to model this situation. Here we implement this model for 
a time-dependent magnetic field, augmented by the sudden approximation to model the 
rapid field ramps. 

By performing a frequency rescaling, an initial condensate density can be achieved 
for two bodies that is comparable to the density range studied experimentally. The 
premise of our model is that the energy scale of the trap energy levels is very low in 
present day experiments, far smaller than the molecular binding energies of interest. 
At the same time, we anticipate that the physics of any single molecule formation 
is controlled by the interaction of just two atoms, even in a many-atom condensate. 
Accordingly, we consider just two atoms in an oscillator trap of very high frequency, 
adjusted so that the density of the two atoms becomes the same as the condensate 
density in the experiments. The resulting approach is then used to model the recent 
experimental results of Donley et al. pT. Our two-body model is shown to describe 
most of the nontrivial features observed in the experiment, although some discrepancies 
remain. This may indicate that a more elaborate inclusion of the many-body effects 
may be necessary to achieve a full quantitative description. Nevertheless, our results 
show that more of the key effects can be interpreted in terms of two-body physics than 
appears to have been realized in existing theoretical models. 

This paper is organized as follows: Section 2 determines the eigenstates of two 
trapped atoms using a quantum-defect-style method that differs from the treatment of 
Ref. [15 but is equivalent. Section 3 discusses the behavior of the atom pair close to a 
Feshbach resonance, and our approximate solution of the time-dependent Schrodinger 
equation using the sudden approximation. Section 4 discusses the frequency rescaling 
employed to interpret the many-atom system. Section 5 compares the results of our 
model with recent experiments, while Section 6 summarizes our conclusions. 

2. Two Interacting Atoms in a Trap 

We consider two atoms of mass m in a spherical oscillator trap of angular frequency u, 
which interact through a zero-range potential ^(r) |18j . 



Here Osc is the two-body atom-atom scattering length and r is the relative coordinate of 
the two particles. The Hamiltonian of the two-body system separates into a center of 
mass part and a relative part. The center of mass part has the usual harmonic oscillator 
solutions which are not affected by the scattering length; hence, we focus on the relative 
motion in the following. 




(1) 
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Since the contact potential in Eq. (Q) acts solely on the s-wave symmetry, we 
consider only solutions of the relative Schrodinger equation with zero orbital angular 



momentum. We define the harmonic oscillator length Oho = yJi/(um/2) corresponding 
to the reduced mass, m/2, as our length scale, leading to a dimensionless radial 
coordinate x = r/oho- Our energy unit is chosen to be hu, resulting in a dimensionless 
energy variable e = E/huj. The s-wave eigenfunction (f)e,i=o of the radial Schrodinger 
equation in the relative coordinate corresponding to energy ehu is then rescaled, 

u{x) 1 



€.1=0 



(2) 



so that the radial equation has only second derivatives, 

(-^^ + ^x^)^x) = e«(x). (3) 

The contact potential, Eq. ((H), imposes a boundary condition on the logarithmic 
derivative of u{x) at the origin: 

u{0) asc ■ ^ ^ 

Solutions of these equations have been obtained by Busch et al. [T^. However, we now 
rederive these solutions in a slightly different manner, along the lines of quantum defect 
theory (QDT) 

We start with a pair of solutions of Eq. 0, / and g, that have regular, 

Ux) = A,xe-^'/\F,{-u-^-x'), (5) 

and irregular, 

g^(x) = B^e-^'/\F,{-u-^;^;x'), (6) 

behavior at the origin, at any energy. Here u denotes a quantum number, u = e/2 — 3/4, 
while Al, and are constants that will be determined later. In the following, our 
solutions are characterized by the subscript u instead of the subscript e. We calculate 
the asymptotic behavior of the two solutions using the known behavior of the confluent 
hypergeometric function, iFi, 

as well as the gamma reflection formula: 

sin(7ri^) 



mT{l-u) = -^. (8) 



For X oo, we obtain 

^ A^T(-) f-x-2--2e-V2sin(7,^,)£(^^±l) + x^'^+ig-"'/^ cos(7rz/) , ^ , 
2 \ n r(z/ + 



2' 
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Our goal is to recast these solutions in the "usual" QDT form [12] given by 

U ^ -C (D-ie"'/2a;-2.-2 gi^^^^) _ jj e-xV2^2.+i cos(7rz/)) , (9) 

g,^C (D-^e"'/2a;-2--2 cos(7rz/) - D e-"'/V"+^ sin(7rz/)) . (10) 

In addition, we want to normalize the functions / and g such that their Wronskian 
^[fuydu] is 2/tt. These requirements can be fulfilled by defining A^, and B^, [Eqs. (0) 
and Q, respectively] appropriately, which leads to 



9.^ 



I^^^^±4xe-^'/\F,{-u;-;x'), (11) 

r Z/ + 1 1 IV , 2' ^ ' ^ ^ 



Using these normalizations, the constants C and D [Eqs. Q and ()10j)] become 
C = 1/0F and D = 0F/-^r(z/ + l)r(z/ + 3/2). 

With these solutions for / and g in hand, we can proceed in the spirit of QDT by 
deriving a solution to the radial Schrodinger equation that accounts for an additional 
non-oscillator short-range potential. For distances beyond those where the short-range 
potential is non- negligible, the radial wave function Ui, must assume the form 

Uu = fu cosTTyU — g^, sinvTyU. (13) 

Armed with the known asymptotic behavior of / and g, we can determine the asymptotic 
behavior of u and impose the requirement that u is finite at large x (i.e., the coefficient 
of the growing exponential is zero). This leads to the equality sin7r(z/ + /x) = 0, which 
can be recast quantization condition 

e = 2(n-^) + ^, (14) 

where n = 0,1, 2,.... The last step towards finding the energy levels e of our confined 
atom pair is to impose the boundary condition on the solution u, which is implied by 
the contact potential at the origin [Eq. (j3))]. Using the fact that at small arguments 
iFi(a;6;x) approaches 1 [201, we obtain 

u'jO) _ f'jO) cos7r/i _ _aho ^^1^^-^ 

u(0) (7(0)sin7r/i 

Upon inserting the explicit forms for / and g as given in Eqs. (jll|) and (|12p. we obtain 
the following equation for the quantum defect fi (see also [T7] ) 

tan7r/i = 7^ . (16) 

«ho r (t + i) 

Equations (jl4j) and (jl6|) allow determination of the energy spectrum for any value of 
the scattering length age- It can be shown that these equations are equivalent to the 
transcendental equation 
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of Busch et al. [14j. The energy quantization conditions derived here for a zero-range 
pseudo potential using QDT also apply to a confined atom pair interacting through an 
arbitrary (i.e., non-contact) short-range potential, the only difference being that the 
quantum defect has a different value. Regardless of the specific short-range potential, 
the corresponding eigenfunctions Uy are then given outside the potential range in terms 
of the hypergeometric U function pTl I14j . 



where Ny is a normalization constant. 

3. Two atoms near a Feshbach Resonance. Overlap matrix elements 

To apply our formalism derived above to two-atom states that lie energetically near a 
Feshbach resonance, we rewrite our scattering length function of the magnetic 

field strength B 



Here, Obg denotes the scattering length far from the resonance, Bq denotes the resonance 
position, while A is a parameter related to the width of the resonance. In our 
numerical calculations (see Sec. E)), we choose parameter values in agreement with 
those obtained by fitting to recent experimental data |22I of the Feshbach resonance 
in ^^Rb: Bq = 155.041 G, A = 10.71 G, and a^g = -443 a.u. Note that a similar 
model was employed in Ref. JH] for the case of a Feshbach resonance in ^^Na. While 
that study additionally incorporated an energy dependence of the scattering length, the 
energy dependence of asc(-B) in the present study is comparatively weak over the range 
of energies and magnetic fields considered. Hence we choose to neglect it. 

Our goal is to describe our two-atom system when it is subjected to time-dependent 
magnetic field ramps B{t), like the ones used in the experiments of pp. To approximate 
the experimental B{t), we assume a piecewise constant B{t) as indicated in Fig. IHl In 
our approach the assumption of an instantaneous variation of B{t) amounts to the use of 
the sudden approximation. In this approximation, the wavefunction of the atom pair is 
initially unaffected by the instantaneous change of B{t). Let the magnetic field strength 
before (after) the instantaneous change be Bi {B2). It is then convenient to write the 
time-dependent superposition state at field strength Bi in terms of the eigenfunctions 
{(pui}ui=Q,oo (with corresponding eigenvalues {eui}ui=o,oo), and that at field strength B2 
in terms of the eigenfunctions {4>y2}u2=o,oo (with corresponding eigenvalues {eu2}u2=o,oo), 



The "new" expansion coefficients a^^ can be expressed through the "old" expansion 
coefficients. 





vl/(S„t)=^a(^)(t)0,^,, J = 1,2. 



(18) 



a 




(1) 
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where 0,^2,1^1 denotes the overlap matrix between the eigenfunctions corresponding to 
Bi and B2, respectively, 



'1^2,1^1 



V2 I Yi^i 



) 



(19) 



The following derivation shows how the eigenstate transformation projections 
0^2^1,1 can be determined analytically. Let Mj^^ j be the reduced radial wavefunctions 
corresponding to 



V\,2 

If we multiply the equation for Uy^ by u^,^ and that for u^^ by 
subtract the resulting equations and integrate the result from to 00, an integration 



by parts gives 



{x)dx 



(20) 



Here, W [uy^, u^^^^q denotes the Wronskian of u^^ and u^^ evaluated at the origin. The 
above formula allows the determination of both the normalization constant [Eq. Q] 
and the overlap matrix element 0^^^^-^ [Eq. (HH)], 

3 

2 cos(7r z/i) r(- + ui) r(l + z/2) sin(7r z/2) 



O 



1^2,1^1 



— 2 cos(7r z/2) r(l 



z/i) r(- + z/2) sin(7rz/i; 



/[4 7r (z/2-z/i)] 



and 

Ny = 



ra + z/' 



r(^ + ^) 



2n + 



sin (2 



TT Z/ 



-1/2 



(21) 



(22) 



where ip denotes the digamma function [20j. The overlap matrix elements Oy2,ui are 
expressed in terms of the z/ quantum numbers, which, in turn, are directly linked to the 
energy eigenvalues e through e = 2z/ + 3/2. Figure Q shows examples of the dependence 
of the overlap matrix elements on the non- integer quantum numbers z/i and z/2. While 
the overlap matrix elements oscillate rapidly with Az/ (Az/ = z/2 — z/i), our analytical 
formula for the matrix elements results in stable numerical calculations (see Sec. Ej). 



4. Description of many-body effects through a frequency rescaling 

The simple model described in the previous section allows us to describe the states of a 
trapped atom pair that undergoes sudden changes of the interatomic interaction, here 
parameterized accurately through a magnetic field- dependent scattering length. Our 
goal is now to apply our two-atom model to interpret an ensemble of atoms (A^ ^ 2). 
Without constructing a rigorous many-body approach like, e.g., the ones based on the 
field theory formalism, we will attempt to account for many-body effects by using our 
two-particle model with a rescaled frequency. The idea suggested by Cornell |2S1 is to 
capitalize on the importance of the diluteness or gas parameter nasc^, {n is the density 
of atoms) which in more rigorous many-body models of degenerate Bose gases plays a 
key role in determining the behavior of the system. Instead of modeling an iV-atom 
system we will model a two-body system that has the same diluteness parameter as the 
experimentally studied atom sample. 
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Figure 1. Dependence of the overlap matrix element 0^2,1^1 between two 
eigenfunctions corresponding to quantum numbers Ui and 1/2 = vi + Ai/ on the quantum 
number difference, A;y, for three values of 1^1: 1^1 = —5 (solid line), = 1.5 (dashed 
line), and Vi = 10 (dashed-dotted line). Note that the solid line in this figure represents 
the projections of the molecular state onto the other trapped atom states. 



To motivate that our two-body description can, at least to a crude level of 
approximation, account for many-body physics we calculate the overlap integral for 
one particular case. We consider the overlap integral between a two-body state located 
very far from the resonance centered at Bq, Uh.o{r), with that corresponding to a value 

of B close to Bq, UM{r). Who refers to the trap ground state, which we approximate 
through a state describing two independent atoms with zero scattering length, 

^ho(^) ~ 4i^c--'/2 (23) 

Um{i^) denotes a molecular state (i.e., the state that remains bound even in the absence 
of the confining potential). Neglecting the influence of the confining potential we assume 
the wavefunction of this state to be: 

UM{r) ^ V2^c-^\ (24) 

where k = l/flsc- We estimate the overlap integral of these two states, by assuming that 
the exponential in Mho is approximately 1 over the range relevant for the evluation of 
the integral. This aproximation yields 




The absolute value of this overlap matrix element gives the probability p that the initially 
"unbound" atom pair ends up in the molecular state as the B field is tuned close to 
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P=4=(-] (26) 



resonance, 



For small p {p <^ V^)) can extend our two-body treatment to model an ensemble 
of N atoms. As the B field is tuned to resonance, each atom pair in the ensemble 
has the probability p to form a molecular bound state. Using this simple picture, the 
calculation of the fraction of atoms that transform into molecules, fatoms^moiecuies, due 
to the magnetic field ramp amounts to a simple counting of atom pairs, 

_ 2 N{N - 1) ^ 

Jatoms^molcculcs ,7 „ p ~ i\ p — . \^ ' ) 

N 2 y/TT aho-^ 

Within our estimate, the fraction of atoms transformed to molecules is proportional to 
the diluteness parameter nasc^ (the density is proportional to N/aho^). This back-of- 
the-envelope estimate, although not a rigorous proof, provides a somewhat quantitative 
motivation for the idea that lies behind the frequency rescaling introduced above. 

Since in our two body model we will choose a scattering length that is the same as 
in the sample used in the experiment, having the same diluteness parameter for the two 
amounts to having the same density. This condition, that our two atoms confined by 
the harmonic potential have the same peak density as the average many-body density, 
n, provides the criterion for determining our rescaled frequency u'. For the purpose 
of evaluating u' from this condition we neglect the atom-atom potential and consider 
the two non-interacting atoms in the ground state of the harmonic trap. Based on this 
assumption we find the relation between the density of the many-body sample, n, and 
the trap frequency confining the two atoms, u', to be 

^ = ( 2 ) • (28) 

''^atom V-^/ 

The experiment jTj uses an inhomogeneous trap that has a geometric mean of ^ 2 tt 12 
Hz, and involves approximately = 17100 atoms. The starting value of the magnetic 
field of about 162.2 G corresponds to a scattering length of 220 a.u.. For these 
parameters, using either the Thomas-Fermi formula or the Gross-Pitaevskii equation, 
one obtains a density profile with an average of approximately n ~ 4 • 10^^ cm~^ (see, 
e.g.. Fig. 10 of Ref. 0). Note that this is lower than the claimed experimental density 
of n = 1.1 ■ 10^'^ cm~^. We find that the final results (see Fig. |3)) depend sensitively on 
density, and this way of rescaling the frequency in Eq. (j^Hjl seems to best reproduce the 
experimental results. The rescaled frequency that we obtain for this density in = 4-10^^ 
cm-3) is u' ^ 2 7r3.72 kHz. 

Figure 121 shows the eigenvalues of two ^^Rb atoms confined by a harmonic potential 
with angular frequency u' as a function of the magnetic field strength B in the vicinity 
of the Feshbach resonance described through the parameters given in Sec. El Notice 
that the sequence of avoided crossings of energy levels associated with the Feshbach 
resonance becomes smeared out at this high frequency. We group the two-body states 
into three groups (see Fig. IS}: the molecular state, the trap ground state, and the 
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155 160 165 

B(G) 

Figure 2. Energy eigenvalues for two ^^Rb atoms for a rescaled frequency uj' = 
2 7r3.72 kHz (corresponding to = 17100) near a Feshbach resonance centered at 
Bo = 155.041 G with Obg = -443 a.u. and A = 10.71 G. We distinguish three groups 
of states: the molecular state (solid line), the trap ground state (dotted line), and the 
trap excited states (dashed lines). 

excited states of the trap. To connect our two-body study to the A^-body systems 
studied experimentally, we make the following correspondence between the two-atom 
states and the many body states. The molecular state of two atoms is in the ground 
state of the center-of-mass motion, and accordingly we associate this population with 
translationally cool (condensed) molecules in the A^-atom system. At fields above the 
resonance, the lowest positive energy state (which we refer to as the trap ground state) 
corresponds to condensate atoms. Finally, the higher trap excited states correspond 
to non-condensed atoms, i.e., the experimentally observed "jets of hot atoms" p. In 
the next section, we interpret the occupation probability of the two-atom states as the 
fraction of atoms ending up in the corresponding many body states. 

5. Numerical Results 

We now use the above description to simulate the experiment of Donley et al. pp. 
This experiment consists of applying two magnetic field pulses, separated by a time 
interval of variable length (denoted by tevoive in Fig. EJ, to a condensed sample of ^^Rb 
atoms. The two pulses of durations ti and t2, respectively, ramp the magnetic field to 
Bm, a value close to the Feshbach resonance, which couples the atomic and molecular 
states of the sample. After the apphcation of these magnetic pulses, three distinct 
components are identified experimentally: a remaining sample of condensed atoms, a 



S 

w"" 



A Two-Atom Picture of Coherent Atom- Molecule Quantum Beats 



11 



hot burst of atoms, and a "missing" component which is not detected under the current 
experimental conditions. At least two theoretical approaches jHl E] ? both involving a 
field theory formalism, identified this third component (i.e., the missing component) 
as a molecular condensate, however, other interpretations exist. The size of the three 
components oscillates as a function of the evolve time tevoive with a frequency that 
corresponds to the energy of a weakly-bound, vibrationally excited state of the ^^Rb2 
dimer. 

We model the experimental situation using the sudden approximation to describe 
the sharp rises and drops of the magnetic field. At time t = in Fig. El our initial state 
is chosen to be the trap ground state corresponding to the initial value of the magnetic 
field, Bi. This initial state is then propagated in time. During the time period where 
the magnetic field is unchanged, B{t) = Bi, the time propagation simply modulates the 
phase of the initial state. After the magnetic field changes to the value B{t) = Bm, 
the system is projected onto the new eigenstates at that field. We then propagate the 
quantum amplitudes of the two-atom energy eigenstates. At the end of the double 
pulse, when the magnetic field becomes Bj, our final state is expanded in the (p^, states 
corresponding to the final value of the magnetic field, Bf. Using the correspondence 
between two-body and many-body states discussed in the previous section, we then 
relate the final population probabilities to the Donley et al. experiment. In taking 
the absolute square of this final quantum amplitude for each distinguishable final state, 
cross terms arise that exhibit quantum beats, the most prominent of which is between 
the molecular state and the atomic trap ground state (condensate). The parameters 
entering our simulations are the rescaled frequency cu', the four magnetic field values Bi, 
Bevoive, Bm and Bf, as well as the three time periods ti, tevoive and t2- Our particular 
simulation values are given in the caption of Fig. El 

Figure |3] shows the results of our time-dependent two-atom calculation for a scaled 
frequency u' that gives a two-atom density of 4 ■ 10^^ cm"^. This figure can be directly 
compared with Fig. 6 of Ref. |IJ. Inspection of these two figures shows reasonable 
agreement in the average values and the oscillation amplitudes for the three components. 
In addition, our results reproduce approximately the relative phases of the oscillations of 
the three different populations. The observed agreement is remarkable if one considers 
that we have used the sudden approximation for magnetic field ramps that in the 
experimental conditions last about 15 /is. We have also assumed that the atomic sample 
has a constant (average) density across the atom cloud, while an improved model could 
consider averaging over a distribution of densities (see, for example, [H]). The frequency 
of the beats that we observe is approximately 150 kHz in accordance with the bound 
state energy produced by our parametrization of the scattering length given in Eq. 
(121). In contrast, the experiment measures oscillations at approximately 196 kHz for the 
same value of i?evoive- In fact, a possible way to improve our model would be to use 
a better description of the atom-atom scattering properties, and consequently, a more 
accurate value of the molecular binding energy. Our estimations show that inclusion of 
the energy dependence up to the effective range level (the first non-zero correction, see 



A Two-Atom Picture of Coherent Atom- Molecule Quantum Beats 



12 



166 



164 



162 



o 



160 



158 



156 



154 




Figure 3. Magnetic field B as function of time t in our modeling of the experiments of 
Donley et ai. The thick dashed line represents the position of the Feshbach resonance, 
-Bo = 155.041G. The abrupt variations of the magnetic field reflect our use of the 
sudden approximation. The parameteres of the pulse are: Bi = 162.2 G, Bm = 155.5 



G, Be 



159.85 G, Bf = 162 G, and ti = 12 ^s, ^2 = 13.6 fis. 




10 15 20 25 30 35 40 



''evolve ^1^^) 

Figure 4. Population of the ground trap state (filled circles) and of the excited trap 
states (open circles) at the end of the magnetic field pulse shown in Fig. 0] The sum 
of the two populations is plotted with squares and allows determination of the number 
of atoms that made the transition to molecules by subtraction from the total number 
of atoms {N ^ 17100). 
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Figure 5. Evolution of population of the three types of states as a function of 
time during one magnetic field double-pulse (with iovoivo = 20 fis, identical to the one 
plotted in Fig. Sohd line: trap ground state population (condensate); dashed line: 
molecular state population; dotted line: population of the excited trap states ("burst" 
atoms). 

[T7] ) significantly improves the agreement between the bound state energy produced by 
our simple model and by a multichannel calculation. However, in order to maintain our 
simple description of the system dynamics during the pulse evolution, we do not include 
this correction here. 

Note that our rescaled frequency yields hoo'/ks ~ 180 nK, which implies that the 
first excited pair above the trap ground state would share approximately two times this 
energy. We note that the resulting energy per atom for the least energetic of our "hot 
atoms" is approximately 180 nK above the trap ground state. This is not very far from 
the experimentaly observed energy (150 nK) of the burst atoms. However, including 
all excited trap states, our model predicts that burst atoms will have a higher average 
energy than in the experiment. In addition, our model predicts a dependence of the 
energy of the hot atoms on the initial atomic density like the one given by Eq. ()28p 
while the experiment does not mention an observable density dependence of the energy 
of the hot bursts. We conclude that the rescaling used in our model allows a reasonable 
prediction of the number of the hot atoms, but has a limited ability to account for their 
energy distribution. 

Figure El shows the evolution of the populations during only one pulse (with 
^evolve — 20 /xs) as a function of time. According to our calculations, the molecular 
eigenstate and the excited states (i.e., the non-condensed atoms) acquire roughly 
comparable populations after the complete first pulse (after t = 22 /is) and mantain 
them until the last ramp of the pulse, when the population of the excited states becomes 
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considerably larger. The significant population we observe for the molecular state 
during tgvoive supports our interpretation of the quantum beats seen in the end-of-pulse 
populations. These beats are the result of the interference of quantum paths that go 
through the intermediate molecular state with those that go through the intermediate 
ground or excited trap states. (The difference between the latter is too small to 
show up on the time scales considered here.) The predictions of this two-atom model, 
concerning the intermediate time populations of the molecular and hot components, are 
consistent with the comments of Braaten et ai, fO] (see also [21]) which point out the 
difficulty of interpreting the intermediate populations shown in Ref. 0. Interpretation 
of the results of Kokkelmans and Holland 8] is complicated by the fact that their two- 
body representation does not consist of eigenf unctions of the molecular Hamiltonian. 
Accordingly, their "molecular state" is not actually the two-body molecular eigenstate, 
except at magnetic fields well above 160 G. In problems with a linear Schrodinger 
equation, a simple basis change could always be carried out, to re-express the physics 
in an eigenrepresentation. Here, however, the nonlinearity of the coupled equations 
complicates this transformation. At the same time, as Braaten et al. fU] comment, the 
success of the final calculations in reproducing the experimental observations with no 
adjustable parameters is immediately apparent and convincing that the right two-body 
physics has been incorporated into the formulation. We mainly recommend caution 
in interpreting the meaning of the "molecular state" in the Kokkelmans and Holland 
formulation, except at high magnetic fields where it approximately coincides with a 
two-body eigenstate. 

The two-atom model disagrees with a specific qualitative prediction of Mackie et 
al. [Hj. Whereas in our model the "hot" atoms at the end of the pulse can be created 
through any of the three intermediate states (i.e, molecules, condensed atoms or non- 
condensed atoms), in the model of Ref. JT] the "hot" atoms are solely the results of 
the "rogue dissociation" of intermediate molecular states. Our approach suggests that 
around half of the final hot atoms are not produced by rogue dissociation. In fact, 
the pathway condensate —>■ hot atoms, which is apparently neglected by Ref. is 
of comparable importance. This may in fact be the additional loss mechanism that is 
cited as being "missing" in Ref. [TTj . 

6. Conclusions 

We have investigated the results given by a simple model of the dynamics of two trapped 
atoms near a Feshbach resonance. Our model accounts for the interaction between the 
two atoms using a zero range potential and it also includes the confinement of the 
atoms by an external harmonic trap. This model can be used to make predictions 
regarding atomic Bose condensate driven with the help of magnetic fields near a 
Feshbach resonance, if a rescaled frequency is used to achieve a density comparable to the 
density of the many-atom condensates studied experimentally. A correspondence can be 
drawn between the three components observed in the recent field-ramp experiments and 
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three groups of two-body states in our model. The experimentally-observed burst of hot 
atoms appears in our model as atom pairs excited to states energetically higher than the 
trap ground state level. The two-atom model is able to predict the populations of these 
three states, in fair agreement with the experimental observations. Accordingly this may 
provide a useful alternative view of the physics of coherent atom-molecule coupling in a 
condensate. The success of a two-body description may initially seem surprising, because 
the dominant physical processes occurring in condensate experiments near a Feshbach 
resonance are normally viewed as being inherently many-body in nature. Nevertheless, 
the present study suggests that a two-body picture, with minimal modifications, is 
sufficiently realistic to be used for simple estimates at a qualitative or semiquantitative 
level. 
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